home *** CD-ROM | disk | FTP | other *** search
/ Pascal Super Library / Pascal Super Library (CW International)(1997).bin / MATH / NRPAS13 / BESSY.PAS < prev    next >
Pascal/Delphi Source File  |  1991-04-29  |  382b  |  19 lines

  1. FUNCTION bessy(n: integer; x: real): real;
  2. VAR
  3.    by,bym,byp,tox: real;
  4.    j: integer;
  5. BEGIN
  6.    IF (n < 2) THEN BEGIN
  7.       writeln('pause in BESSY - index n less than 2'); readln
  8.    END;
  9.    tox := 2.0/x;
  10.    by := bessy1(x);
  11.    bym := bessy0(x);
  12.    FOR j := 1 TO n-1 DO BEGIN
  13.       byp := j*tox*by-bym;
  14.       bym := by;
  15.       by := byp
  16.    END;
  17.    bessy := by
  18. END;
  19.